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P . ABSTRACT 

A previously undetected {Lx < 10^^ erg s~^) source in the strongly star- 
forming galaxy M83 entered an ultraluminous state between August 2009 and 
December 2010. It was first seen with Chandra on 23 December 2010 at 
Lx ~ 4 X 10^^ erg s~^, and has remained ultraluminous through our most re- 
W ! cent observations in December 2011, with typical flux variation of a factor of 

MH ! two. The spectrum is well fitted by a combination of absorbed power-law and 

"-^ ! disk black-body models. While the relative contributions of the models varies 

I ! with time, we have seen no evidence for a canonical state transition. The lumi- 

Vh ! nosity and spectral properties are consistent with accretion powered by a black 

^ ! hole with Mbh ~ 4O-1OOM0. In July 2011 we found a luminous, blue optical 

counterpart which had not been seen in deep HST observations obtained in Au- 
gust 2009. These optical observations suggest that the donor star is a low-mass 
star undergoing Roche-lobe overflow, and that the blue optical emission seen 
cn ' during the outburst is coming from an irradiated accretion disk. This source 

^ I shows that ultraluminous X-ray sources (ULXs) with low-mass companions are 

ff^ I an important component of the ULX population in star-forming galaxies, and 

O 



provides further evidence that the blue optical counterparts of some ULXs need 
not indicate a young, high-mass companion, but rather that they may indicate 
X-ray reprocessing. 
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Introduction 



Point-like X-ray sources with luminosities exceeding the Eddington limit of normal 
stellar-mass black holes {Lx ^ 3 x 10^^ ergs s~^, assuming isotropic emission) are commonly 
known as Ultraluminous X-ray sources (ULXs). The most luminous class of non-nuclear 
sources in galaxies, they are widely believed to result from some extreme form of accreting 
X-ray binary containing a black hole. However, the nature of both the black holes that power 
the ULXs and the companion stars that fuel them remain enigmatic. Briefly, the principal 
competing models for the primaries are: (a) intermediate (between normal stars and AGN) 
mass black holes (lO^M© < Mbh ^ 10'^ Mq) perh aps formed in the collapsed core of massive 
star clusters — many of which are present in M83 (jColbert fc Mushotzkylll999l ): (b) "normal" 
stellar-mass blac k holes (Mbh ^ 20Mq) that are accreting well above their Eddington limit 



( Begelmanll2002l ) and/or that have beamed emission [e.g 
2006|); and (c) "heavy" stellar black hole s (Mrh ~ 30-70M^) 



King et al.ll200ll : iBegelman et al. 



perhaps formed from di- 



rect collapse of raetal-p oor, massive stars (JBelczynski et al.ll2010l : IZampieri fc Robertsll2009 
PakuU fc Mirionil I2OO2I ) — that are accreting near or just above their Eddington limit. 



A number of lines of evidence suggest that most ULXs are associated with young, 
high-mass stellar populations — the extreme of the high-mass X-ray binary (HMXB) pop- 
ulation. Statistically, ULXs are foun d mostly in star-forming spirals or irregular galaxies 
(llrwin et al.l I2004J : ISwartz et al.l 120041) and t h e number of ULXs per galaxy increases with 
the star-formation rate ( Colbert et al.l 12004 : iLiu et al.l l2006). The cumulative luminosity 
function for ULXs in star-forming galaxies is c onsistent with the extrapolation of that for 
HMXBs JWalton et al.ll201ll : ISwartz et aPboilh . Within individual galaxies, ~75% of ULXs 
are found in thin spiral arms and dust lanes, and tend to have high absorbing column densi- 



^Based on observations made with NASA's Chandra X-ray Observatory, the NASA/ESA Hubble Space 
Telescope, Swift , the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, and the Gemini 
Observatory. NASA's Chandra Observatory is operated by Smithsonian Astrophysical Observatory under 
contract # NAS83060 and the data were obtained through program GOl-12115. The HST observations 
were obtained at the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc. (AURA), under NASA contract NAS 5-26555. The new _ff5T observations 
were obtained through programs # GO-12513 and GO-12683. Data in the HST archive from program GO- 
11360 was also used. The ground-based observations were obtained through NOAO which is operated by 
Association of Universities for Research in Astronomy, Inc. for the National Science Foundation. 
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ties (JLiu et al.ll2006[ ). Theoretical considerations su ggest that high-mass do nors can maintain 
a ULX in a persistent high state for up to 10 Myr (JRappaport et al.ll2005[ ). 



Yet there is mounting evidence for a second population of ULXs, a population that 
represents an extreme form of low-mass X-ray binaries (LMXBs). The cumulative luminosity 
function of ULXs in elliptical galaxies is consistent with th e extrapolation of the luminosity 
function for LMXBs (jSwartz et al.ll2004J : IWalton et al.ll201ll ). suggesting that those ULXs are 
extreme LMXBs. In spiral galaxies, the number of ULXs is not purely a function of the star 
formation rate, but is also a function of the total stellar mass, suggesting tha t 15% to 25% 



of th e total number of ULXs are associated with the older stellar population ( IColbert et al. 
2004J : IlIu et al.l [20061 1 Winter et allbood iMushotzkvl I2OO6I ) . 



Point-like, blue optical count erparts have been identified for a number of ULXs [e.g., 
Ptak et al.ll2006l : iGrise et al.ll2008l ). thus reinforcing the link with a young, high-mass stellar 
population. However, some ULX counterparts that were initially ident ified as massive, early- 
type stars on the basis of their blue colors (e.g. iKuntz et al.l l2005l ) are now thought to 
be low-mass stellar donors and the blue colors due to optical e mission from the accretio n 
disk an d /or repro c essed radiation from the X-ray source (e.g. ICopperwheat et al.M2007l ). 
Indeed, iTao et al.l ( 120111 ) have argued that the optical emission from the majority of ULX 
counterparts is dominated by X-ray reprocessing. However, there have been no unequivocal 
identifications of low-mass donors in such systems, because we have never had a chance to 
observe them in quiescence. 

In this paper we report the dis covery of a ULX that recently erupted in the nearby 
{d = 4.6 ±0.2 Mpc. lSahaet al.ll2006[ l grand-design spiral galaxy M83 (NGC5236), and that 
we can confirm to be powered by a low-mass companion. The first exposures of a large 
Chandra project to study M83, taken on 23 and 25 December 2010, revealed an unexpected 
gift: a ULX — a non-nuclear source with a flux comparable to the total circumnuclear flux — 
that had not appeared in previous Chandra images from 2000-01, nor in any other prior 
X-ray image of M83. The new source appeared ~ 1' east of the nucleus, well away from the 
spiral arms and from the many regions with active star formation. Since its discovery, we 
have monitored the source for a year with Swift, and with Chandra in our ongoing program 
observations. We have also obtained optical images of the field containing the ULX from the 
Gemini South telescope and from the Hubble Space Telescope (HST), and we find that a new 
unresolved source, not seen in previous HST observations of the same field, has appeared 
coincident with the ULX. 

The remainder of this paper is organized as follows: In §2 we present the new Chandra 
and Swift data, along with a brief survey of archival X-ray data where the ULX was absent, 
and in §3 we give the results of spectroscopic and timing analyses. In §4 we describe the 
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optical observations and results. In §5 we argue that the M83 ULX is powered by accretion 
from a low-mass companion that has recently expanded to overflow its Roche lobe, and that 
the present bright optical emission results from the reprocessing of X-rays in the accretion 
disk. We go on to discuss the probable geometry and black hole mass in some detail, and 
close in §6 with a brief summary. 



2. X-ray Observations 

The object was discovered using Chandra, whose spectacular angular resolution allowed 
the quick and definitive determination that no such source appeared in previous X-ray im- 
ages. The source was discovered when its Lo.s-io ~ 4 x 10^^ ergss~\ a ULX by most 
definitions of that class. The Chandra spectrum is unequivocally that of a power-law-like 
source, but there are many classes of sources with similar spectra, so the nature of the 
source was not immediately clear. We turned to Swift to follow the short-term evolution 
of the source, and to archival data from other missions to determine what might have been 
there in the past. 



2.1. Chandra Discovery 

The recent Chandra data were obtained as part of a detailed study of M83 (Program 
12620596; Long P.I.). The data, totalling 729 ks, were obtained in ten observations, each 
longer than 50 ks. The observations are clustered in December 2010, March 2011, and August 
2011, with a final observation in December 2011 (See Table[TJ) We used the back-illuminated 
S3 chip for maximum soft response, since most of the M83 disk fits within its 8' x 8' field. We 
carried out the observations in the "very faint" mode for optimum background subtraction. 
We filtered and analysed the data with standard imaging and spectroscopic tools, su ch as 



dmcopy, dmextract and specextract, in the CIAO Version 4.3 (JFruscione et al.l l2006l ) data 
analysis system. 

As shown in the upper right panel of Figure [H the ULX appears ~ 1' east of the nucleus 
of M83. A count-weighted mean of the centroid positions over the first nine observations 
gives a position 

R.A. (2000) = 13^37°'05M35 ± 0'.014, Dec. (2000) = -29°52'0/'2 ± 0'.'3 

(90% confidence level). The root-mean-square scatter of the source positions derived from 
each Chandra observation is a ~ 0'.'25. We also determined a source position using only the 
five observations during March and April 2011, in which the aimpoint was closest (~ 76") 
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to the ULX, when the point spread function (PSF) would be most narrow and symmetric, 
with a 90% encircled energy radius ~ (/.'Q; we obtained the same result as from the total av- 
erage, with a root-mean-square scatter a ~ O'.'IS. We confirmed the accuracy of the Chandra 
position further by using two sources in the S3 chip with both X-ray and radio detections: 
SN1957D (located ~ 2f3 north-east of the nucleus) and a background radio galaxy [th V 
north-west of the nucleus). The mean X-ray positions of those two sources are within k, 0'.'2 
of the radio positions, measured from our Australia Telescope Compact Array observations 
at 5 Ghz; the radio positions have themselves an uncertainty of ^ 0'.'2. We conclude that 
the mean Chandra coordinates can be offset by no more than Ri (//S from the true position. 

The two earlier observations of M83, taken with the Chandra ACIS-S in 2000 and 2001, 
do not show the pre sence of the ULX (see the upper left panel of Figure [T]). Applying the 



Bayesian method of iKraft et al.l (Il99ll ) to the number of detected counts in the source and 
background regions, we estimate the net count rate in 2000-2001 to be < 1.0 x 10^^ ct s^^ 
at the 90% confidence level. This corresponds to an emitted luminosity ^ 1 x 10'^^ erg s~^ 
for any of the spectral parameter values found in the 2010-2011 series of observations. The 
newly erupted source had brightened by at least a factor of 3000. 



2.2. Swift Monitoring 

After the discovery of the ULX in the Chandra data, we monitored it with several series 
of short (si 3 ks) observations with the Swift X-iaj Telescope (XRT), to follow its short-term 
evolution. 

The ULX is somewhat over 0'.5 outside the extended emission from the M83 bulge, in a 
region free of bright, strongly structured diffuse emission, but in a region where the diffuse 
emission is still significant, as can be seen in the lower right panel of Figure [H Thus, we had 
to choose the background region for photometry from the Swift XRT, which has a pixel size 
of 2'.'4, carefully. We set the source region to have a radius of 18", which contains 70% of the 
encircled energy at 1 keV, and the background region to be an annulus stretching from 18" 
to 35", avoiding several nearby regions of enhanced diffuse emission. The background region 
contains the extended wings of the ULX, comprising some 14% of the total source flujo. 
Examination of the Chandra image for the source and background regions show that both 
contain several faint point sources; none of these is readily detectable in the Swift image. 



^ The XRT instrument handbook provides formulae for calculating the encircled energy at 0.5 keV and 
4.0 keV. We used these formula to calculate the encircled energy fractions for 0.3-2.0 keV (using the values 
for 0.5 keV) and for 2.0-10.0 keV (using the values for 4.0 keV). 
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suggesting that the error produced by not excluding them is on the order of the Poisson 
statistics of the background. We extracted counts from the source and background regions, 
solving the simultaneous equations to get the total source counts and the background counts 
per pixel. Since the exposures were relatively short, the background rate was somewhat 
uncertain, particularly in bands narrower than the full 0.3-10.0 keV range of the detector. 
In order to reduce this uncertainty, we determined the mean background rate over all the 
observations and recalculated the source rates. The difference among individual background 
rates was small, and the signal-to-noise of the source counts was improved by using the time- 
averaged background. The resultant light-curve is shown in Fig. [2l the count rate varies from 
^ 0.03 to ^ 0.07 count s-\ 

The (2.0-10.0 keV)/(0. 3-2.0 keV) hardness ratios, shown in Figure [31 vary significantly 
among the observations, making the use of a single energy-conversion factor inappropriate. 
Instead, we calculated the fluxes and hardness ratios expected as a function of power-law 
index, assuming an absorbed power law with Nh = 1.2 x 10^^ cm~^, determined from the 
best-fitting parameters for the December and March Chandra data. Then, we used the 
measured count rates and hardness ratios to produce a flux for each Swift observation. The 
Swift light curve in flux units is shown in the right-hand panel of Fig. [21 Even with 3 ks 
exposures, the total number of counts is insufficient for spectral fitting, and the variation in 
the hardness ratio argues against summing different observations. We obtained three Swift 
exposures that were nearly simultaneous with Chandra exposures. The fluxes calculated 
from the Swift count rates and hardness ratios are roughly consistent with those derived by 
fitting the Chandra spectra, after the latter were corrected for pile-up. 

The only earlier Swift XRT observation useful for this work was made in 2005. It 
provides an upper limit of Fx(0.3-10 keV) < 1.7 x 10"^^ ergs cm'^ g-^ (Lx(0.3-10 keV) 
< 4.5 X 10^'^ ergss"^), assuming an absorbed power law with N^ = 1.2 x 10^^ cm~^ and 
r = 2.0. 



2.3. Previous X-ray Observations from Other Missions 

M83 has been a popular target in surveys of normal galaxies. As a result there are 
sufficient data to determine whether the ULX in M83 has been bright in the past. 

XMM-Newton: Three observations of M83 are available in the public archive. The 
galaxy is well centered in the field of view of the first of these, from 2003, while in the following 
two, from 2008, the galaxy falls on a peripheral chip, so the ULX is not covered by all of the 
instruments in the last two observations. We measured upper limits to the EPIC count rates 
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using a source region with r < 18" and a background region with 18" < r < 35"; we converted 
count rates to fluxes assuming an absorbed power law with N^ = 1.2 x 10^^ cm"^ and 
r = 2.0. If we assume that the source was constant over all the XMM-Newton observations, 
the upper limit becomes Fx(0.3-10.0 keV) < 3.9 x IQ-^^ ergs cm" ^ s'^ (Lx(0.3-10.0 keV) 
< 1.0 X 10'^^ ergs s~^), while the best individual instrument /exposure result is for the PN in 
2003 of Fx (0.3-10.0 keV) < 7.3xl0-^5ergs cm'^ s"^ (Lx (0.3-10.0 keV) < 1.8xl037ergs s'^). 

ROSAT: There is one PSPC observation from 1993 and two HRI observations from 



1993 and 1994 in the archive, each with exposure time ~ 24 ks ( llmmler et al.l 119991 ). There 
is no obvious source at the location of the ULX in the HRI exposures; the detection limit 
for the combined HRI exposures is Fx(0.1-2.4 keV) si 10~^^ergs cm~^ s~^ (Lx(0.1-2.4 keV) 
K, 2.5 X 10'^^ ergss"^). The nucleus and the ULX are not resolved in the PSPC image. 
However, since in its ultraluminous state it is comparable to the nucleus in brightness, had 
the source been ultraluminous during the PSPC observation, one would expect the central 
source to have had an east-west elongation, which was not observed. Thus we may be 
reasonably confident that the source was not ultraluminous during the PSPC observation. 
The same argument can be applied to the ASCA observation in February 1994, and the 
Einstein IPC observation in July 1979. 

Einstein: Besides the IPC observation from 1979, with an exposure time of 6 ks, there 
are two HRI observati ons from 1980 and 1981 with exposures of 25 and 20 ks, respectively 



(ITrinchieri et al.lll985[). Ta king the dimmest detected source in the co-added HRI images 



from iTrinchieri et al.l (Il985[ ) as an extreme upper limit, we find that the flux is Fx < 1-3 x 



10 ^^ ergs cm ^s ^ {Lx < 3.3 x 10^^ ergs s ^). 

Summary: M83 was observed to a depth that would have revealed a ULX, or even a 
moderately bright {Lx > 10^^ erg s'^) X-ray binary in 1979, 1980, 1981, 1993, 1994, 2000, 
2001, 2003, 2005, and 2008. The continuously high emission for at least 12 months, from 
late December 2010 through late December 2011, suggests that if the source has repeated 
outbursts, they are relatively long. Thus the non-detections over the last three decades 
suggest that this source is either new or has a long period between ultraluminous episodes. 



3. X-ray Results 

The Chandra data are magnificent, providing spectra with > 10^ counts, but have 
sparse temporal coverage. The Swift/XBT data provide a more complete light-curve, but 
can provide no more than a hardness ratio. The analysis of each informs and is informed by 
the other. 



3.1. Temporal Variation 



The Swift count rate light-curve in Figure |2] shows that the count rate varied within a 
range of roughly a factor of two. The data hint at a flux decline since March 2011, which 
is in agreement with the Chandra flux trend. However, given that the source showed flux 
variations of comparable amplitude during January and February, it may be too early to 
conclude that the outburst is near its end. The spectral shape, as tracked by the 2.0-10.0 
keV/0. 3-2.0 keV hardness ratio, is also strongly variable. The count rate versus hardness 
ratio plot in Figure [3] suggests two different regimes of behavior. In epochs 2 to 5, as deflned 
in Table [H the hardness ratio varied strongly while the count rate did not, staying within 
~ 0.03 to 0.045 count s~^. In epochs 6 to 13 the hardness ratio did not vary signiflcantly, 
while the count rate did. After epoch 13 the source returned to the mode seen in the earlier 
epochs. The broad-band count rate and the hardness ratio are not well-correlated, and do 
not follow a characteristic track in the hardness-intensity diagram. 

The broad-band flux light-curve derived from a combination of Swift and Chandra 
observations also fluctuates within a range of roughly a factor of two, between 1 and 
2xl0~^^ ergs cm^^ s^^. We determined the Chandra Huxes from multicomponent spectral 
flts. We derived the Swift fluxes by assuming a power-law spectrum absorbed by a column 
of 1.2 X 10^^ cm~^, constraining the photon index from the hardness ratio, and normalizing 
the flux to obtain the measured 0.3-10.0 keV band count rate. The full range of flux was 
spanned by observations in epochs 2 and 3 (a day apart) and by epochs 1 and 2 (9 days 
apart) and 6 and 7 (11 days apart). There is a hint of greater variability in the observations 
before mid-February (epoch 8) than after; the clump of observations in March (epochs 9-11) 
shows little variation in count rate or flux. Unlike the count rate versus hardness ratio plot, 
the flux versus hardness ratio plot shows little structure, and the conversion from Swift count 
rates to fluxes adds scatter to the data. 

Overall, the source variability during our flrst 12 mont hs of observations is more similar 
to the v ariability seen in UL Xs such as Holmberg II X-1 ( iGrise et al.ll2008l ) and Holmberg 
IX X-1 ( Kaaret fc Fend 20091). than to canonical state transitions, rout inely seen in Galactic 
black hole binaries (JFender et al.l l2004l : iMcClintock fc RemillardI |2006| ) . 



We also investigated the intra-observational variability for each Chandra epoch. We ex- 
tracted background-subtracted lightcurves with the CIAO t ask dmextract, an d analysed them 
using the Icstats, ef search and powspec tasks in FTOOLS (JBlackburnl Il995l ) . We computed 
the power spectrum from the shortest period not biased by the Chandra readout rate (~ 10 
s) to the length of the shortest observation (~ 50 ks). We do not flnd signiflcant features in 
the power spectral density over this 0.1 to (2 x 10^^) Hz band for any observation, nor do we 
see any dips, eclipses or flares. The x^ probability of a constant light-curve is ~ 1 for each 
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epoch; the Kolmogorov-Smirnoff probability ( iKolinogorovl 1 1 94 ll ) that the flux is constant is 
> 7% for all observations except for that of September 2011 (when it is ~ 4 x 10~^). We 
can only place a 3cr upper limit of ~ 45% (~ 16%) to the RMS fractional varia tion for 10 
s (100 s) bins at any epoch. Overall, the short-term variability is unremarkable. iHeil et al. 
(120091 ) has demonstrated that ULXs show a wide variety of short-term variability, including 
many that show no signs of variability. 



3.2. X-ray Spectroscopy 

For each Chandra observation, we extracted source spectra from a circular region with a 
4" radius, and background spectra from an annular region between radii of 6" and 12"; we ver- 
ified that the background region contains no other point sources and only a n egligible amount 



of di ffuse emission. We used the specextract script in CIAO Version 4.3 (JFruscione et al. 



20061 ) to build response and area-corr ected ancillary response files, and we fitted the spectra 
in XSPEC Version 12 (jArnaudI 119961 ). Two fully representative examples of the Chandra 
spectra and the fits we carried out are shown in Figure HI 

We fitted the 0.3-10 keV spectrum from each of the Chandra observations with two 
XSPEC models commonly applied to the study of black-hole X-ray binaries: an absorbed 
disk blackbody (diskbb) plus pow er- law model, and an absorbed thermal Comptonization 
(compti) model (JTitarchukI Il994l ) . The diskbb plus power-law model in ULXs is typically a 
good but purely phenomenological test for the presence of a soft excess below 1 keV. Under 
certain assumptions, the soft excess emission may be attributed to the accretion disk and may 
be used to constrain the inner disk size. The comptt model is a suitable test for mild spectral 
curvature, in particular at the high energy end of the ACTS bandpass; a characteristic class 
of ULX spectra with a downturn around 5 keV is forma lly well fitted with Comptonization 
models at low plasma temperatures, /cTg ?a 1.5 to 2 keV (IStobbart et al.ll2006l : lRobertsll2007l : 
Gladstone et al.ll2009l ). 



20 cm~2 flKalberla et al.ll2005f ) while 



The Galactic absorbing foreground column is 4 x 10 
the column density due to M83 in th is direction is 1.5 x 10^° cm"^ (using the naturally 
weighted map from lWalter et al.ll2008l ). Thus, other than any internal absorption, the emis- 
sion should be only lightly absorbed by intervening components. As a result of these con- 
siderations, we assumed a fixed Galactic foreground absorption with a column density of 
4 X 10^° cm~2, and a variable absorption from within M83 and the system itself. These 
column densities are roughly comparable to the optical extinction (see Section 4.2). 

At the observed ACIS-S count rate ^ 0.2 ct s^^ (Tables [2] and [3]), the source spectra are 



-10- 

affected by pile-upj, meaning that there is a high probabihty that more than one photon falls 
in a single pixel (or adjacent pixels) between frame readouts. In general, one can correct 
for the spectral distortions caused by pile-up either by removing pixels with a high pile- 
up probability from the spectral extraction regions, or by using convolution models during 
spectral fitting. We started with the latter technique, and used the pileup convolution 



model of iDavid ( l200ll ) within XSPEC. This model has two principal parameters which can 
be allowed to vary: the grade migration parameter, a, and the fraction of the events in 
the source extraction region to which the pileup will be applied, psffrac. When we allowed 
both of these parameters to vary, we usually found that there were two separate regions of 
parameter space that fit the Chandra spectra. The two model fits implied very different 
pile-up corrected count rates and different ranges of the a and the psffrac parameters, but 
otherwise produced the same spectral shape and goodness of fit, as determined by x^- To 
break the degeneracy and select the proper range for a and psffrac, we extracted another set 
of source spectra from annuli that excluded the piled-up central pixels. We then constrained 
the a and psffrac parameters so that they would produce a spectral fit to the piled-up data 
that was consistent with the fits to the spectra extracted from the annuli. That this rather 
involved process produces the correct flux is supported by the several nearly simultaneous 
Swift observations for which the pile-up corrected Chandra fluxes are roughly consistent with 
their Swift counterparts. 

We find that a diskbb component is significantly detected in five of the observations 
but not the other five (Table [2]). Although a cool disk-blackbody component produces a 
marginal improvement in the x^ fo^' two of those four observations, the F test shows that 
it is not statistically significant (F-test probability P = 0.654 for the 23 December 2010 
observation, P = 0.203 for 4 September 2011). For the 28 December 2011 observation the 
disk-blackbody produces a significant improvement in the fit (F-test probability P = 0.011) 
although the disk-blackbody component itself is only marginally detected (see Table |2]). 
We see no evidence for a disk-blackbody component in the remaining two observations. 
A soft excess is detected significantly in the 3 April 2011 observation (F-test probability 
P = 1.1 X 10^"^), but at very low temperatures {kT < 0.16 keV) in an energy range where 
the ACIS-S3 sensitivity is low and most of the emitted flux is absorbed. Thus, we can 
neither determine whether that soft excess is a diskbb component, nor reliably constrain its 
characteristic temperature and emission radius. 

The absence of a diskbb component in some observations cannot be due simply to shorter 
exposure times. Had such a component always been as strong as during the March 2011 



http://cxc. harvard. edu/ciao/download/doc/pileup^bc.pdf 
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observations, when it accounted for almost half of the emitted X-ray flux, it would have been 
easily detected at every other epoch. It is implausible that an accretion disk would form, for 
example, in the two days between 23 and 25 December 2010, or vanish between 29 March 
and 3 April 2011. Thus, we suggest that the differing significance of the diskbb component 
in various observations more likely depends on the fraction of the inner-disk photons that 
are upscattered in a rapidly variable Comptonizing medium; it is this upscattering that is 
probably responsible for the power-law component. 

As the observations from March 2011 are those where the diskbb component is most 
prominent, they give us the best direct view of the inner disk, and allow us to use those 
observations to constrain the disk parameters. The characteristic peak temperatures at those 
epochs are ~ 0.3 to 0.4 keV corresponding to a characteristic inner radius rin/(cos0)^/^ ~ 700 
to 1000 km. If Tin corresponds to the innermost stable circular orbit of a face-on Schwarzschild 
black hole, this value would suggest a black hole mass ~ 80 to IOOM0, or higher for a 
fast-spinning black hole. On the other hand, the classical Eddington-luminosity argument 
provides a lower limit of ~ 40Mq for the black hole mass for the same data, Lx(0.3 — 
lO.keV) ~ 5 X 10^^ ergs s~^. These two estimates stand in apparent contrast to one another. 
However, we have no reason to assume that we are seeing disk emission from near the 
innermost stable circular orbit. If the inner part of the disk is covered by an optically-thick 
scattering region, the thermal component should come from a larger radius. This is likely 
to be the case if the mass accretion rate is super-critical, that is above the threshold at 
which a radiatively-efficient, geometrically-thin disk cannot survive. At such accretion rates 
(which may be the defining characteristic of ULXs), the inner region may be shrouded by a 
mass-loaded, optically-thick outflow, and thus the observed temperature of the soft thermal 
component would be the temperature at the photosphere of such outflow, while the fitted 
inner-disk ra dius would be the radius at which the disk begins to be covered or replaced by th e 



outflow fe.g. lKing fc Poundsll2003l:IPoutanen et al.ll2007l : ISoriall2007l : [Gladstone et al.ll2009r ). 



It has been suggested (JGladstone et al.l 120091 ) that a ULX with a lower diskbb temperature 



(/cTin < 0.2 keV) and relatively weak soft thermal component compared to the hard power- 
law component will have a stronger outflow and greater mass loading, and thus a lower 
electron temperature in the wind-dominated, Comptonizing inner region ( "warm corona" at 
kTf. ^ 2 keV), which results in a spectral break at energies ~ 5 keV. 

In our case, the power-law component is dominant or comparable to the diskbb emission 
at all epochs. The strong power-law component suggests that the inner disk may indeed 
remain partly shrouded even in the March 2011 observations, and hence the characteristic 
masses and sizes based on r^^ must be taken as upper limits. Fitting the spectra with 
the comptt model helps us investigate this issue. For all observations, the comptt model 
provides statistically equivalent fits to the diskbb plus power-law model. For all except the 
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23 December 2010 observation, we can place only a lower limit (typically, Tg > 2 keV) to 
the best-fitting electron temperature in the scattering region; this is another way of saying 
that the power-law-like portion of the spectrum extends beyond the upper limit of the ACTS 
energy band, without a sharp downturn at least up to ~ 7 keV. To test for the presence of 
a high-energy spectral break, we fitted the 2.0-10 keV portion of the spectrum with both 
power-law and broken power-law models, after fixing the pile-up model parameters. We find 
(Table Hj) that a broken power-law model does not improve the fit in any of the Chandra 
observations. The lack of a break in the power-la w suggests either that this ULX was not 
in a slim-disk/optically-thick warm-corona state ( JRobertd 120071 : [Gladstone et al.ll2009l ). or 
that its characteristic coronal temperature was higher than typically found in that variety of 
ULXs. We conclude that the accretion rate was only moderately super-critical, perhaps not 
high enough to launch a massive outflow, and that the spectral appearance of this ULX is 
closer to the "very high state" of Galactic B H transients than t o the most extreme examples 
of warm-corona/outflow sources described in [Gladstone et al.l (120091 ) . 



With only a few exceptions, both the powerlaw+ diskbb fits and the comptt fits indicate 
that the source has very little internal absorption. The diskbb fits suggest intrinsic absorbing 
column densities ~ a few times 10^*^ cm~^, somewhat larger than the local H I column, 
1.5 X 10^° cm~^; the comptt fits suggest values consistent with no M83 absorption at all. 

Finally, we searched for emission lines, either from the ULX itself [e.g., Fe K lines) or 
from spatially unresolved. X-ray ionized gas in its surroundings. We combined the spectra 
and their respective response files from all the 2010-2011 Chandra observations to increase 
the signal-to-noise ratio of any possible features. We do not find any significant lines or 
edges. A narrow 6.4 keV emission line could be added to the local best-fitting continuum, 
but can have an equivalent width no greater than ~ 12 eV (90% confidence level). Non- 
detec tion of either broa d or narrow fluorescent Fe lines is not a surprising feature for a ULX 



le.g- 



Goad et al. 



20061). although such lines have been detected in the peculiar case of a 



ULX in M82 ( Strickland &: HeckmanI 120071 : iDewangan et al.ll2006f). In general, Fe lines from 



stellar-mass black holes are weaker than those from AGN (JRoss fc FabianI 120071 ) . because 
iron in the inner region of the accretion disk is more highly ionized, and the blurring of the 
line profile due to Compton scattering is more severe. This is the case in particular when the 
power-law-like X-ray emission from a scattering corona is comparable to or stronger than the 
direct thermal emission from the accretion disk — as is the case for most ULXs above 1 keV. 
If the inne r reg ion of the disk is covered or replaced by an optic ally thick scattering corona 



( Roberts! 120071 ). we do not expect to see any reflection features. ICaballero-Garcia &: Fabian 
(J2010l ) proposed that ULX spectra with a steepening around ~ 5-7 keV are dominated by a 
reflection component with a relativistically broadened Fe line; however, this interpretation 
remains somewhat controversial as it requires ad hoc parameters. In any case, the continuum 
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spectrum of the ULX in M83 does not have such high-energy steepening, either. 

In conclusion, had we not had previous observations of M83, the X-ray properties of this 
source would look like those of any number of ULXs. Its power-law photon i ndex fa 2, an d 



its re l atively low disk temp erature are consistent with average ULX properties (JFeng Sz Soria 
20 111 : iBerghea et al.l 120081 ) . Flux variability by a factor of a few o n timescales of days (and 
longer) is typical of ULXs (JGrise et al.ll2008l : iKaaret fc Feng||2009l ). though statistics on that 
variabilty are sketchy, and the most closely monitored sources are brighter than this one. 
Finally, the lack of a disk-dominated phase of the outburst (kn own as high/soft state in 
Galactic black holes) is also a typical ULX property (ISoriall201ll ). 



4. Optical Counterpart 

Optical counterparts of persistent ULXs are typically faint, usually with a B-type star- 
like appearance at extragalactic dist ances. When significa ntly detected, optical variability is 



only a fraction of a magnitude (e.g., IZampieri et al.ll201l[ ) and its physical interpretation is 



unclear. But in this case, given the proximity of M83 and the dramatic change in the X-ray 
flux, we have a unique opportunity to investigate and understand the corresponding changes 
in the optical and UV fluxes. The change was, as expected, below the detectability of the 
Swift UVOT, but was recorded with high confidence by HST and ground-based observatories. 



4.1. Ground-Based Observations 



As part of a previously planned Gemini program, we were able to image the section 
of M83 containing the ULX from the 8.2-m Gemini-South telescope on 8 April 2011 (UT) 
(Program # 2011A-0436, Winkler PI). We used the Gemini Multi-Object Spectrograph 
(GMOS) in its imaging mode, which has a field 5'.5 square and 0'.'0728 pixels. We binned 
the image 2x2 for an effective pixel size of 0'.'146, appropriate for the seeing of ~ 0'.'7 at the 
time of our observations. Four filters were used: u', g', r', and i', with individual exposures 
of 600, 100, 150, and 200 s, respectively, and with four dithered exposures through each 
filter. We processed the image data using standard IRAFq and Gemini reduction procedures 
for overscan and bias subtraction, fiat-fielding, aligning, and stacking. Flux c alibration was 
done using observations of the spectrophotometric standard LTT 4316 (see iHamuy et al. 



"^ IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the AURA, 
Inc., under cooperative agreement with the National Science Foundation. 
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19921 ) taken imniediately after the M83 images. All observations were at airmass < 1.1. 



A blue stellar object, faint but clearly visible, is present at the position of the ULX in 
the stacked images in the u', g' , and r' bands, and possibly in i' as well (Figure El central 
panel). Even though the source is located in an inter-arm region well within the innermost 
spiral arm of M83, the background from the galaxy is complex enough to make accurate 
photometry of a faint object difficult. 

We had previously taken images of M83 with the IMACS instrument on the 6.5-m 
Magellan I telescope in April 2009 as part of a program to inventory supernova remnants 
and other nebulae in M83 (Program # 2009A-0446, Winkler PI). These images were taken 
through narrow-band filters and so are not directly comparable to the 2011 Gemini images, 
but they nevertheless show fainter stars than the 2011 images because of long exposures and 
exceptional seeing, ^ 0'.'4. No object at the ULX position is visible in these 2009 images 
with an upper limit m^i^Q > 25 mag. Although the bands are far from identical, we can 
use the Magellan images (after PSF matching and scaling) to subtract much of the stellar 
contribution in the Gemini images near the ULX, to produce a background that is far more 
uniform and thus more amenable to aperture photometry. Figure shows both a Gemini 
and a Magellan image of the region, together with a difference image that clearly shows the 
ULX counterpart. From the u' and g' images we subtracted the sum of [OIII] plus a 200 A 
wide continuum band centered at 5200 A; from the r' image we subtracted one in a 150 A 
wide continuum band centered at 6800 A. We used these difference images for conventional 
aperture photometry of the ULX counterpart to obtain the Gemini AB magnitudes given 
in Table [51 The uncertainties quoted include any due to possible residuals from imperfectly 
subtracted stars and diffuse background. 



4.2. Hubble Observations 

Images with better angular resolution and higher sensitivity for the M83 field containing 
the ULX were taken in August 2009 with th e newly installed W FC3 instrument on HST as 



part of the Early Release Science program (JDopita et al.l 120101 ). Within the error circle of 



the ULX position, there are no blue sources; the only objects are a few very faint red stars. 
Following the discovery of the ULX counterpart in the Gemini images, we requested and 
were awarded two orbits of Director's Discretionary HST time for a second-epoch WFC3 
observation, which was carried out on 27 July 2011 (Program 12683, PI: Winkler). The data 
were processed with the standard WFC3 pipeline, including drizzling of the three dithered 
frames taken through each of the four filters. The ULX counterpart is readily apparent in the 
2011 WFC3 images through the F336W (f/), F438W (5), and F555W(y) fihers. Alignment 
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and subtraction of the 2009 / image renders the counterpart easily visible on the 2011 F814W 
(/) image as well. The world coordinate systems associated with the WFC3 images at both 
epochs have been updated to agree with that which we determined for the Magellan images, 
which in turn was based on several hundred astrometric stars. The position of the ULX 
counterpart is 

R.A. (2000.) = 13^37"'05.127^ Decl. (2000.) = -29°52'06.92" . 

We estimate the absolute uncertainty in the coordinates to be, at worst, O'.'l, and the regis- 
tration of all the images relative to one another is better than 10 mas. The difference between 
the optical and X-ray coordinates is 0'.'28, which is well within the combined uncertainty of 
the X-ray and optical positions. WFC3 images from both the 2009 and the 2011 epochs are 
shown for a small region directly around the ULX position in Figure [61 

We have carried out photometry of the ULX counterpart and neighboring stars in the 
WFC3 images using DAOPHOT, consistently taking an aperture with a radius of 2 pixels 
(O'.'OS) and using standard encircled-energy corrections to account for the missing flux. In the 
2009 images, there are a number of very faint red stars quite near the position of the ULX 
counterpart, though it is not clear if any of them is actually the donor star in quiescence. 
After updating the world coordinate systems of both images so they are accurately registered, 
we carried out photometry of the 2009 images using the identical positions and apertures 
as in 2011 to obtain the limits given in Table El While we used standard procedures for 
the July 2011 observations, we used "forced photometry" at the position of the ULX for the 
July 2009 data. We did see an excess at the position of the ULX, and these are the values 
reported in Table |5l However, since there was no specific counterpart at this position, the 
values represent upper limits to the flux from the counterpart. 

As part of an ongoing analysis of the HST data for M83, Hwihyun Kim (private com- 
munication) provided photometry and reddening estimates for stars surrounding the ULX 
position. She determined a mean total extinction of Ay = 0.30 ± 0.34. This value indi- 
cates modest but variable extinction, consistent with the appearance of the HST imagery. 
We adopt this mean value for the ULX in correcting the values in Table El to obtain in- 
trinsic fluxes. From the WFC3 data we find that the counterpart's absolute magnitude is 
Mv = -4.85, and a total optical luminosity of L(3300 - 9000A) ^ 2 x 10^^ ergs s"^ at 
the time of the July 2011 observation. In quiescence, the counterpart was no brighter than 
My f« —2.1, and therefore the optical counterpart has brightened optically by a factor of at 
least 10. 

Comparison between the AB magnitudes obtained on 8 April from Gemini and those 
obtained on 27 July from HST (Table (Sj) indicates no evidence for optical variability over 
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this period. Further monitoring of the optical source would be valuable, either in tight coor- 
dination with X-ray observations, or especially if the X-ray object starts to fade significantly. 



Discussion 



We have discovered a bright new source in the low-extinction inter-arm region of M83. 
Although it is well within the D25 radius, it is still worth considering whether the source 
could be a background AGN. Using a standard definition of the X-ray/optical flux ratio 

log ifx/h) = log fx + 5.5 + mR/2.5 



( iHornschemeier et al.ll200l[ ) where fx is the 0.3-10 keV flux and ttir is the Cousins magni- 
tude, we find for this source, with mpj, ^ 24 and log fx ~ --12, log(/x//R) ~ 3.5. This is 
typical of stellar-mass s ystems, while AGN have typical —1 < log(/x//R) ^ 1 (JBauer et al. 



2004 : iLaird et al.l 120091). Furthermore, AG N at this flux level are relatively rare (only ?« 0.1 
per square degree: ICappelluti et al.l l2009l ) and nearly always have identified optical coun- 
terparts. Finally, we are not aware of any AGN that has been observed to vary in X-rays 
by the factor of at least 3000, as seen for the X-ray source in M83. Thus, it is highly im- 
probable that the source is a background AGN. There are no known historical supernovae 
in this region, so it is unlikely that the sudden X-ray/optical increase is due to a previously 
undetected remnant beginning to interact with the circum-stellar medium. 

Our discovery of a transient ULX in an interarm region of M83 suggests that the 
ULX population is more diverse than often assumed. Most ULXs in nearby galaxies are 
variable sources but have been persistently active throughout the years since their original 
discovery (typically, with ROSAT in the 1990s). Instead, this source was not detected by 
Einstein, ROSAT, XMM-Newton, nor in previous Chandra observations. A flux increase of 
> 3000 between the 2000-2001 and the 2010-2011 Chandra observations is very unusual 
for a ULX, but it resembles typical behavior of Galactic black hole transients. Its current 
bright state has lasted at least twelve months, but likely less than two years (based on the 
faintness of its optical counterpart in the August 2009 HST obse rvations). This is longer 
than most Galactic black hole transients (typically, a few weeks: iMcClintock fc Remillard 
20061 ). but is not unp recedented. For example, the 1996-1997 outburst of GRO J1655— 40 
l asted for 15 months (jSoria et al.ll2000l ). and GRS 1915+105 has remained bright since 1992 
(ICastro-Tiradolboili l. 



The other defining characteristic of this ULX is that it is located far from any star- 
forming region, and it must have a low-mass, evolved donor star (mass < 4Mq), since no 
OB counterparts were detected at the ULX position before the start of the outburst (see 
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5.2). Ilrwin et al.l ( l2003l . 120041 ) argued that ULX candidates in the old stellar populations 



of elliptical galaxies were mostly due to background AGN contamination. In spiral galaxies 
with mixed pop ulations, ULXs with low mass counter parts had been identified only in a 
statistical sense ( ISwartz et al.ll2004l : IWalton et al.l 1201 ll ). However, our secure identification 
of this ULX with its optical counterpart confirms the existence of two different classes of 
ULXs. More specifically, this source is a ULX powered by accretion from a low-mass donor 
in an older environment within a galaxy with active star formation. This suggests that a 
classification of ULXs based on the global properties of their host galaxies is incomplete, as 
it may miss or underestimate a population of older, transient ULXs with short active phases, 
particularly at the lower reaches of the "ultra-luminous" lumin osity range. The p resence of 
older ULXs in star- forming galaxies has long been proposed by iMushotzkyI (120061 ). Further 
evidence for the existence of ULXs with low-mass donors in old stellar populations comes 
from two transient sourc es that are only a fa ctor of two less lu minous than the new ULX in 



M83: one in NGC5128 (ISivakoff et al.l 120081 ) and one in M31 flKaur et al.l l20 111 ). 



5.1. The Nature of the Accretor 



Beginning with IShakura &: SunyaevI (Il973l ). a number of steady-state mechanisms have 
been proposed to allow large super-Eddington accretion rates that lead to much milder 
super-Eddington lurainosit ies. Galactic black-hole binaries however, with a few exceptions 
(jjonker &: Nelemansll2004l ). stubbornly remain below the Eddington limit. While it is ac- 
cepted that ULXs are powered by accreting black holes, the main unsolved issue is whether 
they contain a different (more massive) kind of black hole than typical Galactic sources, or 
are simply in a different accretion state {e.g., at a super-Eddington accretion rate). We found 
that the 0.3-10 keV spectrum of this source is typical of ULXs: dominated by a power-law 
with a photon index ^ 2 (intermediate between the soft and hard state of Galactic black 
holes), with an additional soft thermal component at kT ^ 0.3 keV (cooler than typical 
accretion disks of Galactic black holes). The characteristic inner-disk radius implied by the 
thermal component is ~ 1000 km: this corresponds to a Schwarzschild black hole mass 
?« 100 M© if we are directly seeing the disk all the way to the innermost stable circular orbit, 
or < 100 Mq if the inner disk is hidden by a Comptonizing region. The latter scenario is more 
likely, given the dominant power-law component. In fact, in some Chandra observations the 
thermal component disappears altogether, which suggests that a larger fraction of the disk 
emission is Comptonized at such times. 

A strict application of the Eddington limit requires a black hole mass > 40 Mq; however, 
luminosities up to ~ 3 times Eddington (implying correspondingly lower black hole masses) 



-18- 



can be produced in standard accretion scenarios. For example, analytical solutions of stan- 
dard accretion models show that the true emitted luminosity ~ (1 + logM) for M > 1 
( jPoutanen et al.l 120071 ). Recent radiation-magnetohydrodynamic simu lations confirm that 
the is otropic luminosity can reach ~ L^dd during super-critical accretion (jOhsuga fc Mineshige 
201ll ). In addition, mild geometrical beaming can further increase the apparent luminosity by 
a fact or of two for a standard disk seen face-on and by a factor of 10 for supercritical accretion 
flows ( lOhsuga fc Mineshigdl201ll ). Observationally, several neutroi i star X-ray binaries have 



reached luminosities of a few times Eddington in their flaring state ( jBalucihska- Church et al. 



2010l : iBarnard et al.l 120031 : iHoman et al.l 120071 ). Thus, calculated black hole masses may be 
overestimated by a factor of < 3. 

It is still an unsolved theoretical problem whether black holes can form from standard 
stellar evolution at high metallicity ( iBelczynski et al.l l2010l : iHeger et al.l 120031 ) . as metal 
lines make mass l oss more eff i cient, decreasing the pre-collapse stellar mass drastically. 
As pointed out in iHeger et al.l ( l2003l ). there are very large uncertainties on exactly what 
mass/metallicities can ac t ually produce a black hole . It has been argued (jPakuU fc Mirioni 
2OO2I : iMapelli et al.l l2009l : IZampieri fc RobertsI |2009|) that ULX formation e ither requires a 



low-m etallicity environment, or is enhanced by a low- metallicity environment (jPrestwich et al. 



201 ll ). This source demonstrates that low metallicity need not always be the case, since the 
environment in the inner disk of M83 where the ULX is located can hardly be charac- 
terized as a low-metallicity one. Accurate measurements of abundances in regions with 
relatively high metallicity are notoriously difficult, but recent estimates based on deep 
spectra of H II regions within the D25 radius give oxygen abundances for a radius of IfO 
from the nucleus to be 12-|-log('O/H)=8 .73±0.01 (JBresolin et al.l l2009l ). cor npared to the 



solar value of 8.69 ( jAsplund et al.l l2005l ) and the local ISM value of 8.69 (ISnow fc Witt 



19961 ). A different method applied to the H II regions of M83 produces for the same radius 
12-|-log(O/H)=8.59±0.01 (jPilyugin et al.ll2006l ) compared to their local ISM reference value 
of 8.50. While there is some disagreement about the absolute oxygen abundance, it is clear 
that the relative oxygen abundance for this region is slightly higher than solar. In such an 
environment, even an older population is unlikely to have a low metallicity, unless the ULX 
is from a disrupted dwarf galaxy. This object suggests that black holes with masses > 40 Mq 
can be found in environments where the local abundance is solar to somewhat super-solar, 
wherea s models suggest that at these abundances only black holes with < 15 Mq are pro- 
duced (IBelczynski et al.ll2010l ). Reducing abundances to LMC values can produce black holes 
< 30 Mq, so finding a > 40 Mq black hole in a region with solar abundances, even accounting 
for enrichment since its stellar formation, is rather unusual and difficult to understand. 
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5.2. The Nature of the Optical Counterpart 



The optical colors of the system in outburst are not consistent with a simple blackbody- 
like spectrum; the flux density, as shown in Fig. [71 decreases from the near-UV to the 
visible band, but increases again towards the 1 band. This behavior suggests the presence 
of at least two components: one peaking in the UV (implying a characteristic blackbody 
temperature > 20, 000 K and a characteristic effective radius < 20Rq) and one in the IR 
(implying a characteristic blackbody temperature < 4, 000 K and a characteristic effective 
radius > lOOi?©). There are (at least) two ways to interpret this situation. One possibility 
is that the IR emission is dominated by the large outer disk, and the UV excess is the hot 
spot or the irradiated surface of the donor star. An alternative is that the UV peak is the 
Rayleigh- Jeans tail of the emission from an irradiated disk and the 1-band excess is the 
Wien tail of the emission from a cool stellar component. In the latter case the disk must 
be truncated or at least shaded from irradiation beyond R ~ 2Oi?0. If the disk is shaded, 
there will be no signiflcant emission from beyond R ~ 20R(7) as the temperature of a non- 
irradiated disk would be 2000-3000 K Jshakura fc SunvaevlEoTsI : IPrank. King fc Rainell2002[ ) 
and would drop rapidly with dista nce. At those temper atures, the disk gas is entirely neutral, 
and opacity drops to a minimum (JFerguson et al.ll2005l ). It is still debated whether the outer 



disk is optically thin or thick to the continuum in that regime (e.g. ICannizzo fc Wheeleiill984 



Hynes et al.ll2005t lldan et al.ll2010l ): however, in either case, its emission should be negligible 
compared to that from the irradiated part of the disk and the donor star. The pre-outburst 
optical spectral energy distribution contains only a red component which we suggest is either 
the donor star or an unrelated red star along the line of sight. When we add to that an 
irradiated disk model, with disk parameters constrained by the Chandra observations of 
March 2011, we recover the optical brightnesses of the 2011 HST observations. Thus we 
argue that UV/optical emission is due primarily to the irradiated disk, with some stellar 
contribution in the near-IR. 



We used the XSPEC model diskir (JGierlihski et al.ll2009l . l2008l ) to model the irradiated 
disk emission for the 23 March 2011 Chandra data. This observation was chosen because it 
was the longest, had the best constrained accretion disk parameters, and was made close to 
the Gemini observations, which we can then tie to the i75'T observations. We verifled that the 
best-fltting disk and power-law parameters obtained with diskir ior the X-ray data are consis- 
tent with the corresponding parameters obtained with our previous flts. The advantage of the 
diskir model over the simpler, more traditional diskbb plus power-law model is that it allows 
a more natural extrapolation into the optical band, as the power-law component does not ex- 
tend beyond the peak of the disk emission. We found that the UV/optical colors of the 2011 
ifS'T observations are best-fltted with a disk truncation radius rout ~ H, 200 rjn ~ 1.4 x 10^^ 
cm ^ 20Rq, and a fraction of X-ray photons intercepted and re-radiated by the disk to 
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be ~ 5 X 10 ^, a very plausible value for irradiated disks in binary systems. Reprocessing 
fractions between pa 10~^ and ~ 10~^ are suggested both by theoretical m odelling; 



Vrtileket al.lll990l : Ide Jong et al.lll996l : King et al.lll997l : 



vations of accretion disks in Galactic black holes (e.g., iHynes et al.l 12002 : 



le.g., 
Dubus et al.lll999) and by obser- 



Gierlihski et al. 



20091). 



The red excesses in both the pre- and post-outburst spectral energy distributions (SEDs) 
can be represented as a blackbody with T ^ 3500 K and log L/Lq ^ 3.5, implying a radius 
~ 15Oi?0. This set of parameters is consistent wit h an AGB star of init i al mass ~ 2.5-4Mf7) 
and an age ^ 200-800 Myr (Padova stellar t racks: ISalasnich et aPboOol : iBertelh et allboOS 



20091 : iMarigo et al.ll2008t iGirardi et al.ll2010l ). We do not have enough evidence to determine 
whether the faint red component comes from the ULX donor star, or from an unrelated 
evolved star that is projected in the sky within 0.1" of the X-ray source, or even from a 
positive fluctuation in the background stellar density. There are several other faint red stars 
in the field (Figure |6]). If some of the red emission seen from the optical counterpart when 
the ULX is in outburst stems from an unrelated star, then the true donor star of the ULX 
would be even fainter and less massive than we have calculated, but our interpretation of 
the system remains valid. Even if the red excess is due to a single stellar companion, the 
parameters given above are descriptive rather than prescriptive; binary evolution is likely 
to have produced a companion whose structure is somewhat different than an isolated AGB 
star. 

The upper mass/luminosity limits that we place on the optical counterpart are not 
inconsistent with population studies of the surrounding region. Hwihyun Kim (private co- 
munication) provided photometry, reddening, and age estimates for the bright stars in a 
region around the ULX, based on an analysis of the color-color and color-magnitude di- 
agrams. However, none of those stars is close to the ULX (r > 60 pc); the population 
surrounding the ULX is dominated by an older and fainter population than is accessible by 
color-color and color-magnitude diagrams, even with HST. 

If we had observed the optical counterpart only in the X-ray luminous state, we would 
have classified it as a main-sequence early- type B star or a blue supergiant, in agreement 



with most other ULXs for which a uniqu e optical counterpart has been proposed (ITao et al. 



201 ll ). However, in at least one case ( jFeng fc Kaaretl l2008l : [Roberts et al.l 120081 ) . it was 
already suspected that the optical brightness is strongly affected by X-ray irradiation and 
reprocessing on the stellar surface and/or the outer accretion disk. Here, for the first time 
we can prove that this is the case because this source is a transient and we have two sets of 
observations in which the optical brightness is strongly different. It is plausible that some 
other ULXs with apparently massive donor stars in fact have a low-mass donor and a bright 
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accretion disk. 



5.3. The Geometry of the System 



The high mass-accretion rate required by the X-ray luminosity, coupled with the low 
soft-X-ray absorption observed in the X-ray data, imply mass transfer via Roche lobe overflow 
(as should be the case for most ULXs), and hence the formation of a large accretion disk 
around the black hole. We have argued that most of the optical emission in outburst is 
likely to come from the irradiated disk. The required fraction of intercepted X-ray photons 
is comparable to the typical values inferred for Galactic LMXBs, for which the emission is 
thought to be quasi-isotropic. We would not expect the same level of reprocessing if most 
of the X-ray luminosity was beamed along the axis perpendicular to the disk plane. Thus 
the presence of strong optical reprocessing is further evidence that a ULX is not strongly 
beamed. 

Quantitative estimates of the outer disk size and temperature based on broad-band 
flts to the X-ray/optical spectral energy distribution are somewhat model-dependent; the 
disk-blackbody component of the X-ray emission varies from epoch to epoch. However, we 
can obtain useful order-of-magnitude estimates of the system size by applying the empir- 
ical relation between optical brig htness of an irradiated disk and binary period found by 
van Paradiis fc McClintoc^ Jl994J ): 



Ml 



V 



1.57- 1.51 log P^ - 1.141og(Lx/LEdd) 



where Ph is the binary period in hours. In outburst, the absolute visual brightness of the 
ULX counterpart in the standard (Vegamag) system is My ^ —4.9 mag (Table 0, after 
correcting for line-of-sight extincti on). If the X-ray lu minosity Lx ~ (1^3) x LEdd, as most 
classes of ULX models suggest (e.g. iFeng &: Soriall201ll ). this corresponds to an orbital period 
~ 360 to 830 days. Si nce the donor s tar is flUing its Roche lobe, PhP^^^ ~ 10.5 for mass 



ratios 0.01 ^ Q' ^ 1 (jEggletonl Il983l ). where p is the average stellar density in g cm ^. 
Thus, we obtain 3 x 10~^ < (p/g cm~^) ^,1-5 x 10~^. These are ty pical densities of red 



giants or AGB stars ( ISalasnich et al.l l2000t iBertelli et al.ll2008l 120091 ). consistent with our 
interpretation of the system. As a comparison, our independent estimate of a mass ~ 3Mq 
and radius ~ I5O-R0 for the red star apparently associated with the ULX implies a mean 
density p ~ 10^^ g cm~^, in agreement with the expectations if that is the true donor star. 

The size of the accretion disk remains uncertain. Using our lowest estimate for Mbh ~ 
40 M0 and our upper estimate for the Mgtar ~ 4 Mq, the radius of the primary's Roche 
lobe is at least a few hundred Rq. The maximum size of an accretion disk is ~ 80% of 
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the size of the primary's Roche lobe (jPaczynskil 119771 : IWhitehurstlll988l ). which corresponds 
to a characteristic size ~ 10^'^ cm for this ULX. On the other hand, we have argued that 



the optical colors and luminosity in outburst suggest an irradiated disk size ~ 10^^ cm. 
One possibility is that the disk is truncated at that radius, in which case the whole disk 
would be kept ionized by the X-ray illumination. Alternatively, the outer part of the disk 
may be shaded from irradiation, and therefore much colder, neutral, and not a significant 
contributor to the o ptical/ IR emission; the temperature of a non-irradiated standard disk 

19731 ) at i? ~ 10^^ cm would be only ^ 600 K, for the disk parameters 



(^Shakura fc Sunyaev 



fitted to the X-ray emission. We do not have enough information to choose between the 
two scenarios. Whether the outer disk is mostly ionized or neutral has implications for the 
transient behavior. 



5.4. Transient Behavior 

Most Galactic black holes with a low- mass donor star are transient, while most of 
those with an OB donor are persistent X-ray emitters (albeit with flux and spectral vari- 
ability). We suggest that the same may be true for ULXs . A few other examples of ULX 
transients are t wo sources i r i NG C 1365 (ISoria et al.ll2009n. one in M31 (IKaur et al.ll201ll ) 
two il l MlOl kuntz et ahl boosh. one in NGC 3628 Jstrickland et ahlboOlJ ). and one in 
M82 (JFeng fc Kaaretl 120071 ). The transient behavior in accreting stellar- mass black holes 
has been attribute d to thermal/viscous disk instabilities fi Mineshige fc Wheeleiill989l: lOsaki 



19961: 



Lasota 



20011) and/ or to mass transfer instabilities (iHameury et al. 



1986 



Viallet fc Hameurvll2008h . 



19871 . 



1988 



Mass-transfer instabilities operate when the evolved donor star is nearly but not entirely 
filling its Roche lobe, providing only a low accretion rate and thus producing a low-luminosity 
"quiescent" state. The resulting X-ray photons penetrate a few Thompson optical depths 
below the photosphere of the donor star, causing the convective envelope to expand slowly 
until the star makes full contact with its Roche lobe. This contact dramatically enhances both 
the mass transfer rate through the LI point and the X-ray luminosity of the system, which 
is generally seen as an "outburst" state. However, as a larger accretion disk is formed, it 
eventually shades the LI region from X-ray irradiation, decreasing or stopping mass transfer. 
The parameter space for the mass transfer instability in ULXs is largely unexplored. A 
necessary condition is that the irradiating X-ray flux at the surface of the donor star be 
stronger than its intrinsic flux. This is not the case for ULXs with massive donors, but is 
true for this source and probably for other ULXs with low-mass donors. 



In the disk instability model, the disk follows a limit cycle between a hot (higher viscos- 
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ity) and a cold (lower viscosity) state. The difference between the two states is due to the 
sudden jump i n opacity when the te mperature reaches ~ 6, 500 K and hydrogen becomes 
mostly ionized ( ICannizzo et al.lll988f). The instability can be suppressed if the whole disk is 
kept ionized by X-ray irradiation ( iBurderi et al.lll998l : iKing fc Ritterlll998l : iJaniuk fc Czerny 
201ll ). For the observed luminosity of the M83 ULX at the peak of the outburst, the outer 
disk can be kept ionized by irradiation up to a radius Rirr ~ 10^^(/irr/10~^)^''^ cm, where 
/ii-r is the fraction of X-ray flux reprocessed by the disk. From the inferred size of the bi- 
nary system, the accretion disk can exter id beyond this i nstability bounda ry if the outer 
radius reaches the tidal truncation radius (jPaczynskil 119771 : IWhitehurstlll988l ). On the other 
hand, the observed optical colors suggest that the outer disk of the ULX has a temper- 
ature > 20, 000 K and hence is truncated well inside the tidal radius. Thus, we cannot 
tell whether the disk instability operates on the ULX in M83 and will be responsible for 
ultimately bringing the outburst to a close. 



If the outburst does end soon, and the X-ray luminosity declines to below 10 erg s~ , 
we will have a great opportunity to monitor its spectral behavior in the luminosity range 
typical of Galactic black holes. Determining whether or not a ULX behaves like an ordinary 
stellar-mass black hole (e.g., similar state transitions and evolution in the hardness-intensity 
diagram) when its luminosity drops will tell us whether it is powered by an intrinsically 
different type of black hole, or by an ordinary stellar-mass black hole at an extremely high 
accretion rate, not usually reached by Galactic black holes. 



6. Summary 

We have discovered a new ULX in M83 using Chandra, and have characterized the X- 
ray properties of the source in a series of Chandra and Swift observations extending through 
December 2011. We have also detected the optical counterpart to the source using the 
Gemini South Telescope and HST. The ULX is located in an inter-arm region and well away 
from sites of active star formation. Its observed properties and our interpretation of them 
can be summarized as follows: 

1. At its discovery in December 2010, the luminosity was Lx(0.3 — lOkeV) ^ 4 x 
The source has remained bright; its luminosity has varied by a fac- 



]^g39 gj,gg g 1 



tor of two, and it has only recently dropped to ~ 2 x 10 ergs s~ . There is no previous 
evidence for the existence of this X-ray source in observations extending back to 1979; 
the X-ray flux has increased by at least a factor of 3000. 



2. Although there is significant variation between observations, the X-ray light-curves of 
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individual observations show no signs of short-term variabihty, nor are there signs of 
orbital modulations or eclipses. 

3. The X-ray spectra can be well fitted by a disk blackbody plus power-law model, or by 
an absorbed Comptonized spectrum, typical of most ULXs. The strength of the disk 
blackbody component varies on timescales of days. We attribute this to the fraction 
of inner disk photons upscattered in a variable Comptonizing region, rather than to a 
change in the disk. 

4. In those X-ray spectra where a disk is evident, the disk luminosity and temperature sug- 
gest an inner accretion disk radius of about 1000 km, corresponding to a Schwarzschild 
black hole with mass Mbh ~ 100 M0. The black hole mass could be less if the inner 
disk is hidden by a Comptonizing region. In order to strictly obey the Eddington limit, 
the black hole must have a mass Mbh ^ 40 Mq. However, if we allow the possibility of 
accretion luminosities up to three times the Eddington limit at super-critical accretion 
rates, the X-ray data are still consistent with an ordinary stellar-mass black hole. 

5. A blue optical counterpart with My ~ —4.85 has appeared since August 2009, presum- 
ably at the same time as the ULX. The only stars near this site prior to the appearance 
of the ULX are faint and red, so must belong to an older population. The donor star 
is not an OB star, and is likely to be a red giant or AGB star with M < 4Mq and 
age > 500 Myr. We note, however, that had the system been observed only in its lu- 
minous state, the donor could easily have been mistakenly interpreted as an OB star, 
consistent with most other ULXs with uniquely identified optical counterparts. Some 
of these other ULXs may well have low-mass donor stars as well. 

6. During the X-ray outburst, the spectral energy distribution of the optical counterpart 
is dominated by a blue component, with My ~ —4.85 mag, which we interpret as the 
Rayleigh- Jeans tail of the emission from the outer disk, heated by the X-ray photons. 
In addition, there is a faint red component that may arise from the surface of the donor 
star, although we cannot exclude the possibility that it stems from other unrelated stars 
in the vicinity. 

7. The M83 ULX system provides clear evidence that not all ULXs involve an OB donor 
and a young stellar population, confirming the suggestion from statistical studies that 
there are two classes of ULXs. 

8. The existence of a ULX in the inner disk of M83 suggests that it is possible to produce 
black holes also in systems with at least solar metallicity. 



-25- 



The ULX in M83 has provided significant insight into the diversity of the ULX pop- 
ulation, and provides us with an unequivocal example of a ULX with a low-mass donor. 
Continued monitoring of this source will bring greater depth of understanding to this source 
and this class of sources. Determining the length of the current outburst will constrain 
the extent to which sources like this contribute to the persistent ULX population, and the 
length of the decline may constrain the size of the accretion disk. It is still actively debated 
whether ULXs are a state reached by ordinary stellar-mass black holes at extremely high 
(super-Eddington) accretion rates, or are powered by a different (more massive) type of ac- 
creting black hole. Almost all Galactic black hole transients pass through a high/soft state 
(dominated by thermal disk emission) after the peak of their outburst and before returning 
to quiescence; they remain in that state for several weeks. If this ULX does the same, the 
luminosity at which the state changes will help us understand and quantify the relation 
between stellar-mass black holes and ULXs. If, instead, this ULX continues to behave like 
a ULX even when its luminosity goes below ^ 10^^ erg s~^, then it is likely that there are 
intrinsic physical differences between the BHs in ULXs and the ordinary stellar-mass black 
holes. Clearly, this is an interesting object that will continue to illuminate our ignorance 
about ULXs. 
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Table 1. M83 X-ray Observations 



Epoch 


Obsid 


Instrument 


Date 


Exposure 


Flux'' 


Lx 










(s) (10-^* erg cm-2 g-i) 


(10^7 crgs-l) 




Swift 




0005605001 


XRT 


2005-01-24 


8592 


<1.75 


<4.47 


2 


0031905002 


XRT 


2011-01-03 


399 


114 ±33 


290 ± 83 


3 


0031905003 


XRT 


2011-01-04 


1620 


96 ±15 


245 ± 39 


4 


0031905004 


XRT 


2011-01-07 


2213 


192 ± 25 


489 ± 64 


5 


0031905005 


XRT 


2011-01-11 


2140 


139 ± 19 


355 ± 48 


6 


0031905006 


XRT 


2011-01-23 


2896 


183 ± 17 


465 ± 44 


7 


0031905007 


XRT 


2011-02-04 


2938 


113 ±13 


287 ± 35 


8 


0031905008 


XRT 


2011-02-16 


2882 


222 ± 18 


566 ± 46 


9 


0031905009 


XRT 


2011-02-28 


2863 


163 ± 15 


415 ± 39 


10'^ 


0031905010 


XRT 


2011-03-15 


2285 


187 ± 19 


476 ± 48 


11<= 


0031905011 


XRT 


2011-03-24 


3258 


201 ± 16 


513 ± 40 


12 


0031905012 


XRT 


2011-06-25 


3240 


159 ± 15 


404 ± 39 


13 


0031905013 


XRT 


2011-06-30 


3146 


147 ±15 


373 ± 38 


14d 


0031905014 


XRT 


2011-07-27 


3588 


94 ±11 


239 ± 27 


15 


0031905015 


XRT 


2011-08-24 


951 


155 ± 33 


395 ± 84 


16<= 


0031905016 


XRT 


2011-08-29 


2706 


106 ± 15 


269 ± 39 


17 


0031905017 


XRT 


2001-09-04 


4048 


118 ±12 


300 ± 31 


Chandra 




793 


ACIS-S 


2000-04-29 


50981 


< 0.037 


<0.1 




2064 


ACIS-S 


2001-09-04 


9842 






lA 


12995 


ACIS-S 


2010-12-23 


59291 


120-110 


360-300 


IB 


13202 


ACIS-S 


2010-12-25 


98780 


130-120 


440-320 


lOA'' 


12993 


ACIS-S 


2011-03-15 


49398 


150-150 


450-530 


lOB 


13241 


ACIS-S 


2011-03-18 


78963 


160-150 


530-410 


llA'^ 


12994 


ACIS-S 


2011-03-23 


150058 


160-150 


510-420 


IIB 


12996 


ACIS-S 


2011-03-29 


53044 


150-150 


530-410 


lie 


13248 


ACIS-S 


2011-04-03 


54329 


150-150 


500-400 


16<= 


14332 


ACIS-S 


2011-08-29 


52381 


100-100 


290-290 


18 


12992 


ACIS-S 


2011-09-05 


66286 


100-100 


320-270 


19 


14342 


ACIS-S 


2011-12-28 


67103 


80-80 


230-210 


XMM-Newton 




0110910201 


MOSl 


2003-01-27 


21168^ 


< 0.895 


< 2.28 






MOS2 




21670 


< 0.879 


< 2.23 






PN 




9668 


< 0.728 


< 1.85 




0503230101 


MOSl 


2008-01-16 


15946 


not in FOV 








MOS2 




16463 


< 1.90 


<4.82 






PN 




9490 


not in FOV 






0552080101 


MOSl 


2008-08-16 


25413 


< 1.50 


< 3.82 
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Table 1 — Continued 



Epoch 


Obsid 


Instrument 


Date Exposure 

(s) 


Flux'' 
(10-14 ergcm-2 s"!) 


Lx 
(lO^'' ergs-l) 






MOS2 
PN 


26619 
15706 


< 1.65 
< 0.988 


<4.19 
< 2.51 


ROSAT 



rh600024a01 
rh600024n00 
rp600188a02 



HRI 1994-07-30 24156 

HRI 1993-01-20 23529 

PSPC 1993-01-11 23350 



< 1.1 



< 2.54 



ASCA 



61016000 



SIS 1994-02-12 



40186 



Einstein 



I1334S29.XIA 

H1334S29.XIA 

H1334S29.XIB 



IPC 1979-07-31 5708 

HRI 1980-01-15 20002 

HRI 1981-02-13 24574 



< 13' 



< 33 



''In 0.3-10.0 keV. 

''The Swift exposure covered the first part of the Chandra exposure; the Swift exposure lasted from 12:17:32 
to 22:05:56 while the Chandra exposure began at 12:21:40. 

•^The Swift exposure covered the end of the Chandra exposure; the Swift exposure lasted from 11:20:06 to 
21:05:57 while the Chandra exposure ended at 22:18:33. 

^00:34:45 to 03:54:36 This exposure was coincident with the HST exposures in F336E, F438W, F814W. The 
HST exposure in F555W occured shortly after the end of the Swift exposure. 

"The Swift and Chandra exposures were not quite coincident; the Swift exposure lasted from 11:35:00 to 
15:13:56 while the Chandra exposure began at 18:41:51. 

*Aftcr soft proton flare cleaning. 

^The limiting flux is taken from the 3cr limit given by llmmler et al.l 1 119991 ) who comb ined the two HR I 
observations for their analysis. The luminosity has been calculated from their flux using the lSaha et al.l 1120061 ) 
distance. 

''The angular resolution was insufficient to separate the ULX from the nuclear emission. However, it should 
be noted that while in the ULX state, our source has a flux comparable to that of the nucleus, and thus would 
produce a roughly east-west elongation of the nuclear source. This elongation is not seen for this observation, 
suggesting that the source was not in an ultraluminous state. 

'The limiting flux is taken from the analysis of lTrinchieri et al.l l ll985r who combined the two HRI observations 
before analysis. We have used the smallest detected flux fro m a point-like sou rce in that study for our limit. 
The luminosity has been calculated from their flux using the ISaha et al.l 1 1200611 distance. Visual inspection of 
the images reveal no source at the location of the ULX in either image. 



Table 2. Best-fit power-law+diskbb Model Parameters'' 



Parameter 










Val. 


lie 














12995 


13202 


12993 


13241 


12994 


12996 


13248 


14332 


12992 


14342 




Date 


23 Dec 


25 Dec 


15 Mar 


18 Mar 


23 Mar 


29 Mar 


03 Apr 


29 Aug 


04 Sep 


28 Dec 




PSFfrac 


omtlil 


0.581°:;;, 


U.»0_Q Q2 


o.83«;«l 


S4+0-01 
U.»4_ooi 


n 70+0-02 


n 71+0.05 


71+0-15 
'^•'1-0.30 


n 44+0-26 
^•44_o.44 


0A7tr,l 




Ny_{10'^° cm-2) 


4.3li:^ 


9.2tti 


1 l+*-2 
^■1-0.6 


5.4tt-J 


3.4+O-f 


6 1+5-8 
"■1-3.6 


i4.8±t:I 


3.4+2-5 


c: q+4.5 


5.2tr.? 




r 


1 QS+016 

i.y»_g 22 


1.89l«:ii 


, 04+0.39 
^•^4_o.36 


9 ic+O.lS 


9 nq+o-23 

^•Uy-0.23 


2.18±«i| 


9 9^+0.19 
^■^"-0.14 


1.82±«,;» 


1 7fi+0-21 
l-'O-0.16 


i-60tg;ll 




-^Vpl(10"'' phot kcV-i s-i)b 


9 c+0.6 
^■3-0.4 


2.5t^;! 


-, 7+1.9 
^■'-0.9 


2 7+2-3 

^•'-1.4 


1 q+1-0 
1-^-1.0 


3.ii?:^ 


4 0+0.8 


1 s+O-3 
1-8-0.3 


1 S+0-4 
1-8-0.3 


1 1+0-2 
1-1-0.3 




fcTdbb(kcV)=^ 


n r,o+0.13 

^•^•3-0. 23 


o.i8±ni 


n no+0.03 


n OO+0.05 

u.oo_o 08 


n OO+0.04 

U.oo_o,o2 


n 09+0. 06 
^••^^-0.13 






9^^+0-10 
'J-^0_o 08 


9fi+0-15 
U.^0_o 09 




^dbb 


1 l + -'-l 


23.6±f,:| 


2-3l?:? 


3-9li;? 


5.3±?;? 


0+2.1 
o.o_2 






1 t,+15.6 
1-0-1.4 


1 9+16.8 
1-^-1.1 




xVdof 


1.15 


0.99 


0.94 


1.05 


1.02 


1.15 


1.10 


0.92 


0.94 


0.79 




x' 


250.1 


275.3 


221.2 


296.9 


316.4 


244.6 


232.2 


151.2 


195.8 


148.6 




dof 


217 


278 


236 


284 


310 


212 


211 


164 


209 


187 




C (10-1 cnt s-i)d<= 


1 4'^+0-02 


1 c;n+0-01 

t-OU_Q Ql 


1 q-5+0.02 


2.07t^f, 


9 lo + O.Ol 
^•1-^-0. 01 


9 01+0-01 
^•Ul-O.Ol 


1 qn+O-02 

i-y^-0.02 


qs+O-01 
iJ-y^-O.Oi 


1 11+0.01 

1-11-0.01 


f. 00+0.01 
l-S-^-o.oi 


00 


/ (10-12 gi-g f,jjj-2 s-l)c/ 


1 2+0-1 


1 o+O.l 


1 ^+0-1 


i-6tg;} 


1 fi+0-1 


1 5+"-i 


1 5+''-i 
1-0-0.1 


1 0+0-1 
i-'J-o.i 


1 n+o-i 

1-0-0.1 


n s+o-i 


J^ 


L (1039 erg s-l)<=-'' 


o f;+0.2 


A 4+0-5 
4-4-0.4 


4 ^+0-7 


t; + 1.3 
S.d_o.7 


c -1+0.5 


5 3+1-1 


> 5.0 


9 q+O-4 


3.2«:^ 


9 0+0.3 


1 


/dbb^ 


<3% 


(17±I)% 


{Mttr)% 


(36±i?)% 


(5ltIo)% 


(261}^)% 






(7tl)% 


(9±4)% 





''The complete spectral model fitted was pileup X tbabs Gal X tbabs ^x (power-law+diskbb). All uncertainties are for the 90% confidence interval. 

''At 1 keV. 

'^-'^dbb = ['"in(km)/d(10 kpc)] X cos 6 where r^^ is the apparent inner-disk radius and 9 the viewing angle. 

''Before pile-up correction. 

<=In 0.3-10.0 keV. 

'After pile-up correction. 

^Fraction of unabsorbed X-ray luminosity in the diskbb component. 



Table 3. Best-fit comptt Model Parameters'' 



Parameter 



Value 



12995 



13202 



12993 



13241 



12994 



12996 



13248 



14332 



12992 



14342 



Date 


23 Dec 


25 Dec 


15 Mar 


PSFfrac 


0.59l°:li 


n fi4+0-07 


n 70+0-04 
u- 'y-0.05 


AfH(102° cm-2) 


< 1.6 


< 1.0 


7-2^72 


kTo (keV)'' 


n 1 1 +"02 


n 19+0-01 


<0.12 


kTe (keV)= 


1 q+l-2 

^•y-0.2 


> 2.0 


>6.5 


rd 


7 4+0-9 
'•4-0.8 


<8.3 


1 S+5-3 
^•«-0.5 


Xc° 


4,5+1.2 
4-0-1.9 


3-5±i-' 


2.8t: 



18 Mar 

„ _„+0.03 

'J-' ' -0.03 

< 8.0 

15+0-01 

> 1.8 

h2.6 
-2.7 



5.0^ 



5.6 



+6.0 



23 Mar 
n 70+0-02 

< 1.2 

14+0-01 
'^•i^-o.oi 

> 1.9 

hi. 4 
-0.7 
+0.7 



4.9 
6.9 



29 Mar 
0.73 

< 4.9 

0.14 

> 1.9 

5.0 

4.9 



+0.04 
-0.05 



+0.01 
-0.04 



+2.3 

1.1 

+4.3 



03 Apr 


29 Aug 


04 Sep 


0.751°;°^ 


0.68t«i? 


„ 47+0.20 


4.011:° 


<4.5 


<3.5 


1 1 +0-02 
^-11-0.02 


<0.14 


l'^+0-02 
'^■^•^-0.03 


> 2.6 


> 1.8 


> 1.8 


n 7+3.9 

•^■'-0.3 


r, + 16.4 

■^•^-0.6 


7 6+0-7 
'•"-1.5 


0.3tl-° 


O-K' 


2,7+4.2 



C (10-1 cnt s-i)s^ 1.43 ±0.02 
/ (erg cm-2 s-i)h» 1.1^ 
L (103^5 erg s-i)'^J 3.0^ 



1+0.1 

-0.1 

-1+0.1 

J-o.i 



+0.01 

-0.01 

1 2+"-l 

+0.1 

0.1 



1.50 



3.2"* 



1.93 
1.5 
5.3 



+0.02 
-0.02 
+0.1 
0.1 
+0.5 
-1.2 



2.07 
1.5 
4.1 



0.02 
0.02 
+0.1 
0.1 
1.4 
0.1 



2.13 
1.5 

4.2 



0.01 

0.01 

+0.1 

0.1 

+0.2 

-0.2 



2.01 
1.5 
4.1 



+0.01 
0.01 

+0.1 
0.1 

+0.6 
0.2 



1.90 
1.5 
4.0 



+0.02 
0.02 

+0.1 
0.1 

+0.9 
0.1 



0.98 

i.o: 
2.9: 



-0.01 
-0.1 
-0.1 
-0.4 



1.11 
I.o; 

2.7- 



28 Dec 
0.42 

< 3.5 

0.13 

> 2.2 

5.9 

1.0 



+0.29 
0.42 



+0.02 
0.04 



+ 1.7 

1.6 

+4.2 



xVdof 


1.10 


0.95 


0.97 


1.04 


1.05 


1.15 


1.10 


0.93 


0.93 


0.81 


X^ 


239.2 


263.5 


229.5 


294.5 


325.2 


242.7 


231.3 


151.4 


194.8 


151.1 


dof 


217 


278 


236 


284 


310 


212 


211 


162 


209 


187 



+0.01 
-0.01 


0.83lO;°l 


1+0.1 
-0.1 


0.8l°;l 


+0.2 
-0.2 


2.i±;3i 



00 



''The complete spectral model fitted was pileup X tbabs QaiX tbabs ^x (comptt). All uncertainties arc for tiic 90% confidence interval. 
''Seed piioton temperature. 
"^Plasma temperature. 
''Scattering optical depth. 

''Normalization in lO""* XSPEC units (Titarchuk 1994). 
Uncertainty does not converge as the temperature is unconstrained. 
^Before pile-up correction. 
•'In 0.3-10.0 kcV. 
'After pile-up correction. 
J Unabsorbcd luminosity. 



Table 4. Best-fit Power-Law Over the 2-10 keV Band 



Par; 


uiictcr 










Value 












12995 


13202 


12993 


13241 


12994 


12996 


13248 


14332 


12992 


14342 


r 




2.05+0:1° 


1 97+013 


2.50i»:i° 


^■'^1-0.09 


r, yo + 0.06 


9 4c: + 0.14 
^•^''-0.08 


o op+O.OS 


i-84t»:i^ 


1 oc: + 0.16 

i-°>'-o.io 


1 fii+0-33 
1-"1-0.13 


xi 




1.06 


0.90 


0.85 


1.12 


0.97 


1.13 


1.05 


0.91 


0.84 


0.81 


x' 




114.58 


151.62 


108.22 


195.85 


194.40 


117.19 


109.25 


56.4 


86.6 


66.1 


dof 




108 


168 


128 


175 


201 


104 


104 


62 


103 


82 


Ti 




i-99l°;ll 


1 04+0-10 


NA=- 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


Sbr 


(kcV) 


5.311.? 


t).U_2 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


Ta 




3.14lti^ 


4-56±2.4i 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


xl 




1.06 


0.88 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


x" 




111.92 


146.74 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


dof 




106 


166 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


F test 


P = 0.288 


P = 0.066 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


NA 



^Broken power-law models with the constraint that Fa > Fi do not provide any improvement to the xi value. 
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Table 5. Observations of the Optical Counterpart to the M83 ULX 



Date 


Telescope 


Instrument 


Exposure 




Filter 




Magnitude^ 


A(A)b 


Unabsorbed 










(s) 


Band 


Ac (A) 


AA (A) 




(mag) 


Flux {iyY 


2009 Apr 26 


Magellan I 


IMACS 


7 


X 600 


[OIII] 


5007 


50 








2009 Apr 26 


Magellan I 


IMACS 


7 


X 200 


Green Cont'm 


5316 


161 


msiSO > 25.0 






2009 Apr 26 


Magellan I 


IMACS 


7 


X 200 


Red Cont'm 


6815 


216 








2009 Aug 9 


HST 


WFC3 


3 


X 630 


U (F336W) 


3355 


511 


"13355 =28.1±0.8'^ 


0.48 


3.2 X 10-» 


2009 Aug 9 


HST 


WFC3 


3 


X 640 


B (F438W) 


4326 


618 


"14326 = 27.5 ± 0.6^^ 


0.40 


5.4 X 10"* 


2009 Aug 9 


HST 


WFC3 


3 


X 401 


V {F555W) 


5308 


1562 


"15308 = 26.4±0.3'^ 


0.30 


1.3 X 10"^ 


2009 Aug 9 


HST 


WFC3 


3 


X 401 


/ (F814W) 


8030 


1536 


m8030 = 24.4 ± 0.3^^ 


0.18 


7.7 X 10-^ 


2011 Apr 8 


Gemini-S 


GMOS 


4 


X 600 


u' 


3650 


490 


"13650 =23.4 ±0.3 


0.47 


2.4 X 10-*^ 


2011 Apr 8 


Gemini-S 


GMOS 


4 


X 100 


9' 


4780 


1540 


"14780 = 23.8 ± 0.2 


0.35 


1.6 X W-^ 


2011 Apr 8 


Gemini-S 


GMOS 


4 


X 150 


r' 


6340 


1440 


m6340 = 23.8 ± 0.3 


0.25 


1.4 X 10"'^ 


2011 July 27 


HST 


WFC3 


3 


X 608 


U {F336W) 


3355 


511 


"13355 = 23.18 ±0.08 


0.48 


2.0 X 10-6 


2011 July 27 


HST 


WFC3 


3 


X 330 


B (F438W) 


4326 


618 


"14326 = 23.71 ± 0.06 


0.40 


1.7 X 10-6 


2011 July 27 


HST 


WFC3 


3 


X 245 


V {F555W) 


5308 


1562 


"15308 = 23.73 ± 0.05 


0.30 


1.6 X 10-6 


2011 July 27 


HST 


WFC3 


3 


X 487 


/ (F814W) 


8030 


1536 


"18030 = 23.87 ±0.06 


0.18 


1.7 X 10-6 



''AH magnitudes are observed AB magnitudes at the central wavelength. For WFC3, these use the nominal WFC3 zero points given at 
http://www.stsci.edu/hst/wfc3/phot_zp_lbn. 

''Assuming Ay = 0.30 and using reddening corrections from lRomaniello et al.l 1120021) for the HST filters and from lCardelli et al.l 1119891) at the 
central wavelengths for the Gemini filters. 

'^Values obtained from the 2009 HST observations correspond to the excess flux at the location of the ULX and are upper limits to the flux 
of the pre-outburst donor star. 
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Fig. 1. — Upper Left: X-ray image of the inner region of M83, from the Chandra ob- 
servation in 2000. Color coding: red, 0.3-1.0 keV; green, 1.0-2.0 keV; blue, 2.0-7.0 keV. 
The yellow circle encloses the position of the ULX. Upper Right: X-ray image of the 
same region from our 2010-2011 Chandra observations, showing the newly discovered ULX. 
Color coding is as in the upper left panel. Lower Left: Optical color composite image 
of approximately the same region from the 2009 HST/WFC3 Early Release Science data. 
The diffuse red emission is from Ha and U, V, and I bands are shown as blue, green, and 
red. The 30" x 30" yellow box is shown in greater detail in Figure |5l Lower Right: The 
co-added Swift image in the 0.3-10.0 keV band. The source and background regions used 
for the photometry are shown. Each panel is just over 2' on a side and the physical size of 
each panel is ~ 2.7 kpc. North is up. East is left. 
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Fig. 2. — Top: The Swift 0.3-10.0 keV count rate as a function of time. The open 
boxes are the Swift points labeled by their epoch. The filled diamonds are the Chandra 
powerlaw+diskbb spectral fits converted to Swift count rates. The dates of the Gemini and 
i/5T observations are also marked. Bottom: The flux as a function of time. The open boxes 
are the Swift count rates converted to fluxes as described in the text. The filled diamonds 
are the Chandra fluxes derived from the powerlaw+diskbb spectral fits. 
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Fig. 3.— The 2.0-10.0 keV/0.3-2.0 keV hardness ratio as a function of the 0.3-10.0 keV 
band count rate. The scale on the right is the equivalent power law index, assuming a simple 
absorbed power law with an N(H) of 1.2 x 10^^ cm~^. The points are labeled by their epoch. 
The open boxes are the Swift data points labeled by their epoch. The filled diamonds are 
the Chandra powerlaw+diskbb spectral fits converted to the Swift count rates and hardness 
ratios. 
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Fig. 4. — Left: The 23 March 2011 Chandra/ ACIS spectrum fitted with an absorbed diskbb 
plus powerlaw model, convolved with a pileup model. The data points (in red) are binned to 
a signal-to-noise ratio > 7 for display purposes only. Top panel: The separate contributions 
of the two model components [diskbb in orange, powerlaw in cyan). Bottom panel: The 
X^ residuals for the best-fitting model with pile-up correction (red stepped line), as well as 
without pile-up correction (blue data points). See Table 2 for the best- fitting parameters. 
Right: As in the left panel but for the 29 March 2011 spectrum. The data points (in 
green) are binned to a signal-to-noise ratio > 5 for display purpose only. See Table 2 for the 
best-fitting parameters. Note the change in relative strength of the fit components. 
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Fig. 5. — Images showing a 30" square region centered on the position of the ULX (the 3" 
diameter circle) from Left: Magellan IMACS, ([O III] plus narrow green continuum) from 
April 2009, Center: Gemini-S GMOS, (g) from April 2011, and Right: the difference image 
between the Gemini image and the Magellan one after PSF-matching and scaling. North is 
up and east to the left; the dashed box indicates the 10" field shown in Fig. [HI 
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Fig. 6. — HST WFC3 true-color images of the region around the ULX. In all panels, the 
colors are: blue, for the F438W filter; green, for F555W; and red, for F814W. North is up 
and East to the left. Left: Image from July 2011, during the X-ray luminous state; the ULX 
counterpart appears distinctly blue. The field size is 10" square. This panel can be compared 
to the ground-based images in the previous figure. The dashed box shows the region covered 
by the next two panels. Middle: A 2'.'5 by 2'.'5 detail of the July 2011 image. Right: Image 
from August 2009, at the same scale as the previous panel. The bright blue "star" is not 
present; there are many faint red stars near the ULX position in the 2009 image, but none 
is exactly coincident. 
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Fig. 7. — Model of the multi-wavelength spectral energy distribution (SED) in outburst 
and quiescence. The blue curve at short wavelengths is the 0.3-10 keV best fit diskir model 
for the 23 March 2011 Chandra data, which were chosen for this exercise as they are from 
the longest observation with the strongest constraints on the accretion disk parameters. The 
dark blue data points are the pre-outburst i/5'T measurements from 9 August 2009. These 
are actually upper limits for the flux from any actual counterpart. The red curve is the best 
fit assuming that the upper limits to the light from the pre-outburst star actually correspond 
to the light from that star. The magenta data points are from the 27 July 2011 HST 
measurements; the orange data points are from the 8 April 2011 Gemini measurements; and 
the green dotted line is the modeled emission from an irradiated disk with the addition of the 
pre-outburst blackbody. All data and models have been corrected for absorption/extinction. 



